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where: 

VDVj  is  the  VDV  from  a  single  impact  at  a  particular  speed. 

VDVtotai  is  the  total  VDV  which  is  being  calculated 
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1.  SCOPE. 

This  Test  Operations  Procedure  (TOP)  describes  methods  for  evaluating  the  ride  dynamics  or 
ride  quality  of  ground  vehicles  as  well  as  the  vehicle  occupants’  exposure  to  Whole-Body 
Vibration  (WBV).  Ride  dynamics  pertains  to  the  sensation  or  feel  of  the  passengers  in  the 
environment  of  a  moving  vehicle.  Ride  quality  is  a  measure  of  vehicle  mobility  (linking  vehicle 
speed  and  terrain  roughness)  while  WBV  exposure  is  concerned  with  the  adverse  health  effects 
or  serious  injuries  that  may  occur  as  a  result  of  vibration  exposure.  The  technique  for  collecting 
data  to  be  used  for  either  ride  dynamics  or  WBV  exposure  assessments  is  similar.  The  vehicle 
seats  are  instrumented  with  ride  quality  pads  which  contain  accelerometers  molded  in  a  rubber 
disk  to  provide  measurements  in  three  mutually  perpendicular  axes.  The  instrumented  seats  are 
occupied  by  a  vehicle  crew.  Data  are  acquired  while  the  vehicle  traverses  various  test  courses  at 
pre-detennined  speeds.  The  collected  raw  data  are  processed  according  to  three  analysis 
techniques.  For  the  ride  dynamics  analysis,  the  data  are  processed  by  performing  Fourier 
transfonnations  and  creating  frequency  weighted  Power  Spectral  Density  (PSD)  files  for  the 
individual  test  runs.  The  ride  dynamics  analysis  techniques  are  described  in  Paragraph  4.  For 
WBV  exposure  analysis,  the  data  are  processed  via  specifications  called  out  in  International 
Standards  Organization  (ISO)  263 1-1  ’*  and  ISO  263 1-52.  The  analysis  technique  for  ISO 
2631-1  is  described  in  Paragraph  4.2  and  the  analysis  technique  for  ISO  2631-5  is  described  in 
Paragraph  4.3. 

2.  INSTRUMENTATION. 

2.1  Ride  Quality  Pads. 

Vibration  will  be  sensed  by  accelerometers  mounted  in  ride  quality  pads,  which  will  be  placed  in 
designated  vehicle  seats.  The  ride  quality  pads  are  typically  fabricated  in-house  and  contain 
three  uniaxial  piezoresistive  accelerometers  (vertical,  transverse,  and  longitudinal)  mounted  in  a 
rubber  disk.  The  primary  requirements  of  these  pads  are  that  they  should  not  adversely  affect 
occupant  comfort  and  shall  not  significantly  distort  the  buttock-cushion  load  distribution.  These 
pads  generally  confonn  to  the  suggested  design  of  Society  of  Automotive  Engineers  (SAE) 
J10133,  as  shown  in  Figure  1.  The  semi-rigid  disk  is  fabricated  of  molded  rubber  of 
approximately  80  to  90  durometer  (A  Scale). 

2.2  Other  Measurement  Locations. 

Additional  vibration  (acceleration)  measurements  should  be  made  to  provide  a  “dynamic  map” 
from  the  terrain  surface  to  the  seat  (frame  and  floor  of  the  vehicle  are  also  used)  for  correlation 
with  future  testing  that  may  not  have  seat  pad  instrumentation  (e.g.,  endurance  test  vehicles).  As 
a  minimum,  acceleration  measurements  should  be  made  in  the  vertical  axis  at  one  unsprung  mass 
(suspension)  location  near  the  tire,  the  sprung  mass  corresponding  to  the  previous  location 
(frame),  and  the  base  of  any  instrumented  seat.  All  additional  accelerometers  should  be  DC 
coupled  devices. 


*  Superscript  numbers  correspond  to  Appendix  G,  References. 
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Figure  1 .  Suggested  design  for  ride  quality  pad. 


2.3  Vibration  Measurement  Axes. 

Vibrations  are  to  be  measured  utilizing  the  seat  pad  mounted  accelerometers  in  the  vertical  (z 
axis),  fore/aft  (x  axis),  and  lateral  (y  axis)  directions.  The  vibration  measured  along  these  three 
mutually  perpendicular  axes  pass  through  a  point  on  the  interface  between  the  seated  crew 
member  and  the  seat.  This  orientation  is  shown  in  Figure  2. 


az,  ay.  az  -  ACCELERATION  IN  THE  DIRECTIONS  OF  THE  x.  y.  2.  AXIS 
x-axIs  -  DIRECTION  OF  THE  BACK-TO-CHEST 
y— axis  -  RIGHT -TO-LEFT  SIDE 
2  — axis  =  BUTTOCKS-TO-HEAD 


Figure  2.  Measurement  axes. 
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2.4  Test  Subject  Anthropometry  and  Seat  Adjustment. 

Measurements  will  be  made  at  the  driver’s  seat  and  at  other  crew  positions  as  required  by  the 
vehicle  specification  and/or  the  detailed  test  plan.  The  sex,  height,  and  weight  of  each  occupant 
of  the  instrumented  seats,  as  well  as  the  seat  location  within  the  vehicle,  will  be  recorded.  Test 
subjects  should  preferably  have  a  height  and  weight  between  the  Anny  5th  percentile  female  and 
95th  percentile  male,  and  should  not  be  outside  the  height  and  weight  range  allowable  within  the 
Department  of  Defense  (DOD)  Handbook  (HDBK)-743A4,  except  for  unique  circumstances 
(e.g.,  limited  personnel  resources  with  required  skill).  The  US  Army  Anthropometry 
measurements  are  provided  in  Table  1.  The  seat,  if  adjustable,  will  be  properly  adjusted  per 
manufacturer’s  specification.  Subjects  will  maintain  contact  between  the  seat  pad  and  buttocks 
at  all  times.  Loss  of  contact  during  the  test  will  require  a  retest  for  that  condition.  Other  (non- 
instrumented)  seats  will  be  filled  with  an  appropriate  dummy  load  to  simulate  the  weight  of  a 
crew  member. 


TABLE  1.  US  ARMY  ANTHROPOMETRY 


CATEGORY 

HEIGH 

T 

WEIGHT 

centimeter  (cm) 

inches  (in.) 

kilogram  (kg) 

pound  (lb) 

Male  -  Max 

203 

80 

1 13a 

250a 

Female  -  Min 

147 

58 

41 

91 

Male  -  95% 

187 

73.5 

98 

216 

Female  -  5% 

153 

60.2 

49 

109 

a  Age  dependent  -  max  age  (40+)  shown 


2.5  Accelerometers  and  Signal  Conditioning. 

The  pad  and  ancillary  accelerometers  together  with  their  associated  signal  conditioning  shall  be 
capable  of  measuring  root-mean-square  (rms)  accelerations  in  the  0.1  to  100  Hertz  (Hz) 
bandwidth.  An  on-board  or  telemetry  data  acquisition  system  will  be  used  to  acquire  data  while 
the  vehicle  is  operated  on  the  test  courses.  All  data  will  be  digitized  at  a  minimum  of  400 
samples  per  second  per  channel  after  having  been  low-pass  filtered  at  a  minimum  of  100  Hz. 
Ensure  low  pass  filters  are  adequate  to  prevent  aliasing  of  the  data.  If  a  half-round  obstacle 
requirement  must  be  met,  a  vertical  accelerometer  will  be  installed  at  the  base  of  the  driver’s  seat 
to  comply  with  historical  data.  Other  locations  (e.g.,  on  the  seat)  may  be  instrumented  as 
required  by  the  vehicle  specification  and/or  the  detailed  test  plan.  This  (half-round  obstacle) 
accelerometer  will  be  low-pass  filtered  (post  test)  at  30  Hz  (Paragraph  4.4).  Amplifiers  should 
be  adjusted  to  provide  an  acceleration  resolution  of  approximately  0.01  g.  Analysis  parameters 
should  be  chosen  so  that  frequency  domain  resolution  is  0.2  Hz  or  less. 
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2.6  Acceleration  Units. 

Acceleration  units  for  analysis  will  be  meters/sec  for  ISO-263 1-1/ISO  2631-5  computations, 
ft/sec2  for  absorbed  power  computations  and  “g”  for  other  acceleration  computations. 

2.7  Data  Validation. 

Prior  to  analysis,  data  from  each  transducer  and  each  data  run  will  be  checked  for  stationarity 
and  for  errors  such  as  noise,  amplifier  drift,  clipped  data,  etc.  Procedures  used  for  data  validation 
will  be  presented  in  the  test  report. 

2.8  Length  of  Data  Run. 

The  precision  of  acceleration  spectral  estimates  is  related  to  the  length  of  the  data  run.  Data 
should  generally  be  recorded  for  a  period  of  at  least  one  minute  unless  restricted  by  test  course 
length  and  vehicle  speed. 

2.9  Vehicle. 

The  test  vehicle  will  be  described  by  recording  the  information  in  Appendix  A. 

2.10  Vehicle  Speed. 

Vehicle  speed  will  be  recorded  by  an  auxiliary,  calibrated  speed  sensing  system  to  an  accuracy  of 
+  0.2  kilometers  per  hour  (km/hr)  (0. 1  miles  per  hour  (mph)).  The  test  speed  will  be  held 
constant  (to  the  extent  possible)  at  the  target  speed  for  the  duration  of  the  data  run.  The  speed 
Coefficient  of  Variation  (COV  (standard  deviation/mean))  for  each  data  run  should  be  equal  to  or 
less  than  0.10. 

2.11  Test  Courses  for  Ride  Quality. 

Test  courses  used  for  ride  quality  should  have  a  wave  number  spectrum  (power  spectral  density 
function  of  the  instantaneous  terrain  vertical  profile  in  the  spatial  domain)  with  a  slope  of 
approximately  -2  in  the  log-log  domain.  The  rms  roughness  value  will  be  computed  as  the 
square  root  of  the  integral  of  the  wave  number  spectrum  and  will  be  computed  between  wave 
numbers  that  cover  the  frequency  range  of  interest  (typically  0.5  to  80  Hz).  The  affected  wave 
numbers  are  a  function  of  vehicle  speed  (single  wheel  frequency),  and  integration  band  wave 
numbers  between  approximately  0.016  and  2.0  (wavelengths  of  0.5  feet  to  64  feet)  will  suffice 
for  speeds  up  to  40  km/hr  (25  mph).  The  lower  integration  band  wave  number  should  be  reduced 
to  0.01  (wavelength  of  100  feet)  for  speeds  between  40  and  56  km/hr  (25  and  35  mph). 

2.12  Test  Courses  and  Speeds  for  WBV  Analysis. 

Data  must  be  collected  in  a  manner  that  is  representative  of  the  mission  profile  of  the  vehicle  and 
is  typical  of  the  expected  exposures.  Duration  of  measurements  shall  be  sufficient  to  ensure 
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statistical  precision.  Whole-body  vibration  data  must  be  collected  in  accordance  with  ISO 
2631-1  and  ISO  2631-5.  Refer  to  ISO  2631-1  and  ISO  2631-5  for  complete  data  collection 
guideline  details.  In  order  to  properly  assess  health  effects  from  exposure  to  WBV,  the  vehicle 
should  be  tested  under  its  nonnal  operating  conditions.  Vibration  data  will  be  recorded  on  test 
courses  with  vehicle  speeds  and  vehicle  maneuvering  scenarios  analogous  to  the  Operational 
Mode  Summary/Mission  Profile  (OMS/MP).  For  example,  if  the  vehicle’s  OMS/MP  states  the 
vehicle  will  operate  on  primary,  secondary,  and  cross  country  terrain,  then  vibration  data  must  be 
collected  on  test  courses  that  represent  each  of  the  three  terrains.  If  the  vehicle’s  mission  profile 
calls  for  speeds  on  primary  terrain  up  to  50  mph,  then  the  vehicle  should  be  tested  at  increments 
of  10  mph  up  to  50  mph  on  primary  terrain.  Coordination  with  the  vehicle’s  program  office  to 
detennine  test  courses  and  speeds  that  best  represent  the  vehicles  nonnal  operating  conditions  is 
recommended.  If  program  office  guidance  is  not  forthcoming,  Table  B-l  in  Appendix  B 
contains  the  minimum  requirements  for  speed  and  course  terrains  the  vehicle  needs  to  traverse  to 
complete  the  WBV  data  set. 

2.13  WBV  Analysis  of  Operational  Mode. 

Vehicles  that  produce  WBV  in  environments  other  than  traversing  terrain  are  to  be  tested  under 
conditions  that  replicate  the  operational  environment  in  which  the  WBV  is  produced.  For 
example,  backhoes  with  breaker  attachments  generate  large  amounts  of  vibration  while  the 
vehicle  is  stationary.  A  proper  assessment  will  include  WBV  data  collected  while  the  backhoe’s 
breaker  is  in  operation.  Table  B-2  in  Appendix  B  provides  general  examples  of  Operational 
Mode  Analysis. 

2.14  Test  Conduct. 

a.  Acceleration  and  speed  data  will  be  acquired  as  the  vehicle  is  driven  over  each  test 
course  at  multiple  constant  speeds  as  described  in  the  test  plan  for  a  period  of  at  least  1  minute 
(when  not  limited  by  test  course  length)  for  the  ride  dynamics  courses,  and  for  a  sufficient  period 
of  time  to  capture  the  shock  event  when  negotiating  the  half-round  obstacles.  A  speed  increment 
of  approximately  3  km/hr  (2  mph)  should  be  used  because  vertical  absorbed  power  generally 
increases  as  a  power  law  function  of  speed  for  a  given  rms  roughness  (e.g.,  Absorbed  Power  ~ 
k* Speed5).  Following  each  data  run,  the  crew  will  be  polled  to  ensure  that  it  is  safe  to  proceed  to 
the  next  speed.  Each  instrumented  crew  member's  seat  vertical  acceleration  data  will  be 
processed  at  the  test  site  after  each  data  run  to  provide  an  absorbed  power  value.  The  test  on 
each  course  will  be  considered  as  completed  when  the  driver's  vertical  absorbed  power  equals  or 
exceeds  15  watts,  when  any  crew  member  believes  that  an  increase  in  speed  would  create  a 
hazardous  condition,  or  when  the  course  speed  limit  is  reached  (applicable  only  to  relatively 
smooth  courses).  The  limit  value  of  15  watts  was  chosen  to  satisfy  the  requirements  of  at  least 
one  Anny  vehicle  (a  ride  specification  based  on  “9  to  12  watts”)  and  to  enhance  the  ride-speed 
distribution  for  the  health  hazard  assessment  (Paragraph  5).  Data  runs  that  produced  absorbed 
power  values  from  above  1  watt  to  the  maximum  speed  run  will  be  repeated  to  enhance  data 
confidence. 

b.  To  evaluate  the  shock  input  test  criteria,  the  half-round  obstacles  should  be  traversed 
with  either  wheel  or  track  paths  crossing  the  obstacle  simultaneously.  To  evaluate  vehicle  roll 
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dynamics  and  suspension  characteristics,  additional  data  runs  can  be  performed  with  a  single 
wheel  or  track  path  crossing  the  obstacle  while  the  other  wheel  or  track  path  remains  on  level 
ground. 

3.  VIBRATION  EVALUATION. 

a.  Vibration  evaluation  will  be  performed  using  two  techniques;  an  ISO  technique  which 
deals  with  health  effects  of  exposure  to  vibration,  and  an  absorbed  power  technique  which  is 
used  to  describe  speed  limiting  effects  over  rough  terrain. 

b.  The  principal  factors  that  combine  to  detennine  the  degree  to  which  human  exposure 
to  whole  body  vibration  will  be  acceptable  is  described  in  ISO  263 1  - 1 .  Four  possible  effects  of 
vibration  include  Degraded  Health,  Comfort,  Perception  and  Motion  Sickness.  The  frequency 
ranges  of  these  effects  are: 

(1)  0.5  Hz  to  80  Hz  for  Degraded  Health,  Comfort  and  Perception. 

(2)  0. 1  Hz  to  0.5  Hz  for  Motion  Sickness. 

c.  Unless  required  by  the  test  plan,  only  the  issue  of  degraded  health  will  be  evaluated. 
This  type  of  vibration  is  transmitted  to  the  human  body  as  a  whole  through  the  supporting 
surfaces  of  the  buttocks,  back,  and  feet  of  a  seated  person  in  a  moving  wheeled  or  tracked 
vehicle.  Vibration  is  measured  according  to  the  coordinate  system  originating  at  a  point  from 
which  vibration  enters  the  human  body.  The  coordinate  system  for  the  alignment  of  the  vibration 
transducers  is  shown  in  Figure  2.  The  three  principal  areas  of  contact  for  seated  persons  are:  (1) 
the  supporting  seat  surface,  (2)  the  seat  back,  and  (3)  the  feet,  but  only  the  seat  data  are  of 
interest  for  degraded  health  analysis.  Seat  data  are  also  the  basis  for  the  absorbed  power 
analysis.  Vibration  transmitted  to  the  body  through  a  non-rigid  material,  like  a  seat  cushion,  is 
measured  with  the  transducer  interposed  between  the  person  and  the  principal  contact  areas  of 
the  surface.  This  is  achieved  by  mounting  the  transducers  within  a  ride  quality  pad  as  described 
in  Figure  1 . 

4.  ANALYSIS  TECHNIQUES. 

4.1  Absorbed  Power  Technique. 

a.  Another  t  echnique  f  or  e  valuating  hum  an  r  esponse  t  o  vi  bration  or  r  ide  q  uality  i  s  t  he 
computation  of  absorbed  power.  It  is  a  measure  of  the  rate  at  which  energy  is  absorbed  by  a 
human  subjected  to  ride  vibration.  It  is  accepted  as  a  measure  of  human  tolerance  to  vibration 
for  military  vehicles  negotiating  rough  terrain.  The  absorbed  power  for  a  given  location  and  axis 
is  computed  by  multiplying  the  acceleration  power  spectral  density  spectrum  by  the  appropriate 
transfer  function  and  integrating  the  resultant  spectrum.  A  n  advantage  of  this  approach  is  that 
average  ab  sorbed  p  ower  i  s  a  s  calar  q  uantity  an  d  can  b  e  s  ummed  i  n  co  mplex  m  ulti-degree  o  f 
freedom  systems  to  yield  a  single  value  describing  the  total  average  absorbed  power.  Typically  a 
standard  ride  quality  test  procedure  involves  detennining  the  speed  at  which  the  vertical  average 
absorbed  power  reaches  an  upper  limit  of  6  watts  for  different  types  of  terrain  at  selected  crew 


7 


TOP  01-1-014A 
29  November  2011 


location  seats.  The  terrain  is  usually  characterized  by  the  surface  roughness  reported  as  inches  or 
centimeters  rms.  The  vehicle  speed  as  a  function  of  terrain  roughness  obtained  from  this  process 
is  used  as  one  of  the  mobility  limiting  factors  in  the  North  Atlantic  Treaty  Organization  (NATO) 
Reference  M  obility  M  odel.  M  ilitary  v  ehicle  specifications  frequently  in  elude  r  ide  q  uality 
requirements  based  on  absorbed  power  (e.g.,  the  vehicle  must  be  capable  of  producing  a  ride  of 
6  watts  or  less  at  the  driver’s  seat  in  the  vertical  axis  for  a  given  speed  and  surface  roughness). 


b.  The  absorbed  power  in  each  axis  will  be  calculated  from: 


where: 


n 


P  =  E(Q)A2 

1 


P  =  Absorbed  power,  watts. 


A;  =  rms  acceleration  in  ft/sec2  within  the  ith  spectral  band. 


Q  =  KiK0(FiF4  -  F2F3)  /  (F32  +  W  2F42). 


Wi  =  Frequency,  radians/second. 

F  and  K  values  =  Calculated  from  Appendix  C. 

c.  The  result  is  a  frequency  weighting,  scaling  and  integration  of  the  acceleration  power 
spectral  density.  The  weighting  functions  are  shown  (by  axis)  in  Figure  3.  The  factors  have 
been  normalized  to  a  value  of  1 .0  to  show  frequency  without  the  appropriate  scaling. 


d.  For  each  test  course,  the  absorbed  power  will  be  plotted  as  a  function  of  speed. 
Specification  requirements  are  generally  written  in  terms  of  the  driver’s  vertical  absorbed  power 
(only),  but  the  following  technique  will  be  applied  to  other  locations  and  axes,  if  applicable.  A 
non-linear  interpolation  or  a  power  law  (y  =  axb,  linear  fit  in  a  log-log  domain,)  or  nth  degree 
polynomial  curve  fit  will  be  perfonned  and  the  speed  at  which  an  absorbed  power  of  6  watts  is 
achieved,  will  be  interpolated  from  the  data  set.  Absorbed  power  values  of  less  than  1  watt  or 
greater  than  10  watts  will  not  be  used  in  the  curve  fit  process,  if  possible.  The  6-watt  speed  will 
then  be  plotted  as  a  function  of  terrain  roughness  rms  (a  specific  value  for  each  test  course)  to 
detennine  a  ride  quality  curve  for  the  vehicle. 


8 


TOP  01-1-014A 
29  November  2011 


Absorbed  Power 
Normalized  Weighting  Factors 


0  10  20  30  40  50 

Frequency,  Hz 


Figure  3.  Absorbed  power  normalized  frequency  weighting  factors. 


4.2  ISO  2631-1  Technique. 

a.  The  basic  evaluation  method  described  in  ISO  2631-1  utilizes  the  weighted  rms 
acceleration.  The  weighted  rms  acceleration  is  expressed  in  meters  per  second  squared  (m/s2)  for 
translational  vibration  and  radians  per  second  squared  (rad/s2)  for  rotational  vibration.  The 
weighted  rms  acceleration  is  calculated  in  accordance  with  the  following  equation  or  its 
equivalent  in  the  frequency  domain: 


a 


W 


T 


\alit)dt 


o 


where: 

aw(f)  is  the  weighted  acceleration  (translational  or  rotational)  as  a  function  of  time 

2  2 

(time  history),  in  m/s  or  rad/s",  respectively. 

T  is  the  duration  of  the  measurement,  in  seconds. 
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b.  The  crest  factor  may  be  used  to  investigate  whether  or  not  the  basic  evaluation  method 
is  suitable  for  describing  the  severity  of  the  vibration  in  relation  to  its  effects  on  human  beings. 
The  peak  value  is  determined  over  the  duration  of  the  measurement,  i.e.  the  time  period  T  used 
for  integration  of  the  rms  value.  The  crest  factor  does  not  necessarily  indicate  the  severity  of 
vibration.  For  vibration  with  crest  factors  below  or  equal  to  9,  the  basic  evaluation  method 
outlined  above  is  normally  sufficient.  In  cases  where  the  basic  evaluation  method  may 
underestimate  the  effects  of  vibration  (high  crest  factors,  occasional  shocks,  and  transient 
vibration)  an  alternative  technique  should  be  utilized  such  as  the  running  rms  method  or  the 
fourth  power  vibration  dose  method. 

c.  The  running  rms  method  evaluation  method  takes  into  account  occasional  shocks  and 
transient  v  ibration  byuseof  as  hort  in  tegration  time  c  onstant,  f  or  e  xample  1  second.  T  he 
vibration  magnitude  is  defined  as  a  Maximum  Transient  Vibration  Value  (MTVV),  given  as  the 
maximum  in  time  of  aw(t0),  defined  by: 


a 


w 


M 


where: 

aw(/)  is  the  instantaneous  frequency-weighted  acceleration. 
t  is  the  integration  time  for  running  averaging. 
t  is  the  time  (integration  variable). 
t0  is  the  time  observation  (instantaneous  time). 

d.  The  f  ourth  pow  er  vi  bration  dos  e  m  ethod  i  s  m  ore  s  ensitive  t  han  t  he  ba  sic  e  valuation 
method  by  using  the  fourth  power  instead  of  the  second  power  of  the  acceleration  time  history  as 
the  basis  for  averaging.  The  fourth  power  Vibration  Dose  Value  (VDV)  in  meters  per  second  to 
the  power  1.75  (m/s1'75),  or  in  radians  per  second  to  the  power  1.75  (rad/s1'75),  is  defined  as: 

wr* 


where: 

aw(/)  is  the  instantaneous  frequency-weighted  acceleration. 
T  is  the  duration  of  measurement. 
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The  total  VDV  (VDVtotai)  from  separate  events  can  be  calculated  with  the  4-th  root  of  the  sum  of 
the  VDVi  of  the  individual  events  taken  to  the  fourth  power,  or  expressed  as: 

iX 

VDV*  ^ 


The  total  VDV  from  a  repeating  a  specific  condition  can  be  estimated  by  using  an  individual 
VDV;  value  multiple  times  in  the  above  expression.  F  or  example,  the  total  VDV  value  from  N 
impacts  on  a  half  round  at  a  specific  speed  can  be  estimated  with  the  following  expression: 


VDV 


total 


where: 

VDVi  is  the  VDV  from  a  single  impact  at  a  particular  speed. 

VDVtotai  is  the  total  VDV  which  is  being  calculated 

N  is  the  number  of  impacts  for  which  the  total  VDV  is  being  calculated 

e.  The  f  requency  w  eighting  curves  us  ed  f  or  va  rious  di  rections  of  m  easurement  i  n  t  he 
calculation  of  frequency-weighted  acceleration  for  the  seat  surface  as  outlined  in  ISO  2631-1  are 
provided  in  Appendix  D  and  Figure  4. 
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ISO-2631  (1997)  Weighting  Factors 


Figure  4.  ISO  2631-1  frequency  weighting  factors. 


f.  Daily  vibration  exposure  limits  are  determined  from  the  weighted  acceleration,  aw,  in 
each  axis.  The  assessment  is  made  independently  for  each  axis,  but  the  vector  sum  is  used  if  the 
weighted  accelerations  from  two  or  more  axes  are  comparable.  The  evaluation  is  made  using  the 
frequency  weighting  factors  shown  in  Appendix  D  with  the  multiplying  factors  as  shown  in 
Table  2. 


TABLE  2.  MULTIPLYING  FACTORS 


AXIS 

FACTOR 

X,Y 

1.4 

Z 

1.0 

g.  Exposure  times  for  three  conditions  -  no  documented  health  risks,  caution  zone  for 
health  risks,  and  health  risks  likely  are  calculated  from  the  weighted  acceleration  as  follows: 
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2 

(1)  No  documented  health  risks:  Exposure  time  =  1.5/aw  . 

(2)  Health  risks  likely:  Exposure  time  =  6.0/aw  . 

2  2 

(3)  Caution  zone  for  health  risks:  Exposure  time  >  1.5/aw  and  <  6.0/aw  . 

(4)  For  example:  let  aw  =  2.0  m/s  . 

(a)  Exposure  time  =  0.4  hours  for  no  documented  health  risks. 

(b)  Exposure  time  =  0.4  to  1 .5  hours  for  a  health  risk  caution. 

(c)  Exposure  time  =1.5  hours  for  a  likely  health  risk. 

(d)  Thus,  an  exposure  (per  24  hours)  to  this  level  of  vibration  for  less  than  0.4  hours 
should  produce  no  health  risks,  an  exposure  between  0.4  hours  and  1.5  hours  will  result  in  a 
health  risk  caution,  while  an  exposure  of  1.5  hours  or  greater  will  create  a  likely  health  risk. 

h.  The  health  guidance  caution  zone  is  shown  in  Figure  5. 
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Figure  5.  Health  guidance  caution  zone. 
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4.3  ISO  2631-5  Multiple  Shock  Method. 

a.  The  method  conveyed  in  ISO  263 1  -5  is  an  attempt  to  quantify  whole-body  vibration 
containing  multiple  shocks  in  relation  to  human  health.  T  he  assessment  method  described  in 
ISO  2631-5  presumes  an  individual  is  in  good  physical  condition  and  is  maintaining  and  upright 
and  unsupported  seated  posture.  Specifically  this  method  is  meant  to  focus  on  the  lumbar  spine. 
Different  postures  will  result  in  different  responses  in  the  spine.  The  method  attempts  to  create 
predictions  about  the  number  of  years  for  which  an  individual  can  be  subjected  to  a  sample  time 
history  prior  to  the  onset  of  adverse  health  effects.  H  owever,  in  the  introductions  the  authors 
identify  t  hat  t  his  a  ssessment  m  ethod  a  nd  r  elated  m  odels  ha  ve  not  be  en  e  pidemiologically 
validated. 


b.  The  calculation  of  the  spinal  response  acceleration  dose  involves  the  following  steps: 
calculation  of  t  he  hum  an  r  esponse  of  t  he  s  pine,  counting  t  he  num  ber  a  nd  m  agnitudes  pe  aks, 
calculation  of  an  acceleration  dose  by  application  of  a  dose  model  related  to  the  Palmgren-Miner 
fatigue  theory.  The  spinal  response  to  the  seat  bottom  triaxial  input  is  estimated  with  predictive 
models.  The  spinal  response  in  the  horizontal  plane  is  approximated  by  a  linear  single  degree  of 
freedom  lumped  parameter  model.  The  spinal  response  in  the  vertical  axis  is  approximated  with 
a  non-linear  recurrent  neural  network  model.  The  model  coefficients  used  for  the  vertical  model 
are  specific  to  a  s  ampling  rate  of  1 60  samples  per  second.  T he  data  collected  for  this  test  were 
re-sampled  to  160  samples  per  second. 


as: 


c. 


The  acceleration  dose,  Dk,  in  meters  per  second  squared  in  the  k -direction  is  defined 


Dk 


where: 

Aik  is  the  ith  peak  of  the  response  acceleration  aik(t). 

d.  The  average  daily  dose,  Dkd,  in  meters  per  second  squared,  can  be  computed  with  the 
following: 


Dud  ~ 


j— 1  mj 


-i  1/6 


where: 

Dkj  is  the  acceleration  dose  for  the  jth  unique  period  within  a  day. 
tdj  is  the  duration  of  the  daily  exposure  to  condition  j. 
tmj  is  the  period  over  which  Dkj  was  measured. 
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e.  As  s  tated,  t  he  ISO  263  1-5  m  ethod  a  pplies  t  he  P  almgren-Miner  a  pproach  f  or  t  he 
evaluation  o  f  t  he  ef  fects  of  i  ntemal  p  ressure  c  hanges  on  t  he  s  pine.  The  m  ethod  relies  upon 
experimental  data  which  has  shown  the  P almgren-Miner  exponent  varies  for  different  tissues  (5 
to  14  f  or  bone  s,  20  f  or  c  artilage).  F  or  t  he  pu  rpose  of  e  stimating  a  dverse  h  ealth  effects  a 
conservative  e  xponent  of  6  w  as  us  ed.  From  t  his  a  ssumption  t  he  c  onstants  a  nd  s  implifying 
equations  are  developed  to  predict  the  life  cycle  response  of  human  tissue.  The  daily  equivalent 
static  compression  dose,  Sed,  in  megapascals  is  obtained  with  the  Dkd  value  calculated  earlier: 


1/6 


ZKaJ6 


where: 

mx  =  0.015  MPa/(m/s2). 
my  =  0.035  MPa/(m/s2). 
mz  =  0.032  MPa/(m/s2). 

f.  A  factor  R  is  defined  to  account  for  the  repeated  exposure  to  the  same  environment  for 
a  number  of  years.  T  he  method  attempts  to  take  into  account  the  differences  in  initial  tissue 
strength  related  to  the  age  of  the  seat  occupant.  As  the  exposure  time  increases  the  muscle  tissue 
reduces  in  strength.  For  this  reason  R  must  be  calculated  sequentially: 


R  = 


1/6 


where: 

N  is  the  number  of  exposure  days  per  year, 
i  is  the  incremented  year  count  from  the  beginning  of  exposure, 
n  is  the  number  of  years  of  exposure. 

c  is  a  constant  representing  the  static  stress  due  to  gravitational  force, 
b  is  the  age  at  which  the  exposure  starts. 

Sui  is  the  ultimate  strength  of  the  lumbar  spine  for  a  person  of  age  (b+i)  years. 
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g.  The  value  Su;  is  indented  to  account  for  the  variation  in  vertebrae  bone  density,  which 
normally  reduces  with  age.  In-vitro  studies  were  used  to  derive  the  relationship  between  Sui  (in 
megapascals)  and  b+i  (years): 


Sui  =6.75-0.066(6  +  /) 

h.  The  method  provides  some  guidance  on  the  interpretation  of  the  Seci  and  R  values.  For 
instance,  if  the  seat  occupant  is  subjected  to  240  days  of  exposure  per  year  for  a  lifetime,  a  Sed 
value  of  0.5  MPa  indicates  a  low  probability  of  an  adverse  health  effect.  A  Seti  value  of  0.8  MPa 
indicates  a  high  probability  of  an  adverse  health  effect.  These  would  correspond  to  R  values  of 
0.8  and  1 .2  respectively.  Obviously  the  severity  of  Sed  will  depend  on  the  presumed  duration  of 
the  daily  exposure,  tdj,  which  the  shorter  sampling  time,  tmj,  is  intended  to  represent.  Often  the 
daily  exposure  time  is  unknown  and  will  depend  on  the  severity  of  the  duty  to  which  the 
occupant  is  exposed.  Obviously  future  work  would  benefit  from  a  better  understanding  of  real 
world  exposures.  Similarly  assumptions  must  be  made  about  the  starting  age  and  number  of 
years  of  exposure  to  which  an  occupant  is  subjected.  For  an  enlisted  soldier,  an  average  starting 
age  of  2 1  years  with  2 1  years  of  service  can  be  a  good  presumption. 

4.4  Half-Round  Obstacle  Technique. 

The  peak  acceleration  from  the  vertical  accelerometer  at  the  base  of  the  driver's  seat  will  be  low- 
pass  filtered  at  30  Hz  (4-pole  Butterworth  filter  forward  and  backward  to  preserve  phase)  and  the 
resulting  peak  a cceleration  w ill  be  plotted  as  a  f  unction  of  speed  for  each  half-round  obstacle. 
The  s  peed  a  t  w  hich  a  pe  ak  va  lue  of  2.5  g’s  i  s  r  eached  w  ill  b  e  d  etermined  b  y  a  n  on-linear 
interpolation  or  curve  fitting  technique  as  described  above. 

5.  HEALTH  HAZARD  ASSESSMENT  AND  DATA  FILE  STRUCTURE. 


a.  The  collected  vibration  data  should  be  stored  in  simple  text  files.  The  data  collected 
must  b  e  in  B  ritish  C  olumbia  R  esearch  ( BCR)  f  ormat,  w  hich  i  s  a  text  f  ile  w  ith  th  e  e  xtension 
“.BCR.”  The  f  ilenames  m  ust  ha  ve  t  he  “.BCR”  ex  tension  because  it  is  the  onl  y  e  xtension 
recognized  by  the  J olt  software.  The  files  should  contain  two  sections,  Header  and  Data,  and 
should  be  formatted  as  shown  in  Appendix  E. 

b.  The  he  alth  ha  zard  a  ssessment  i  s  ba  sed  on  t  he  vehicle’s  m  ission  pr  ofile,  which  i  s 
reflected  i  n  the  t  est  course  m  ileage  br  eakdown  of  t  he  e  ndurance  t  est.  T  he  s  ame  r  ide  da  ta 
previously  described  will  be  collected  on  the  endurance  test  courses  following  the  guidelines  of 
Paragraph  2.11.  To  improve  test  efficiency,  data  from  remote  test  courses  need  not  be  acquired 
if  the  ride  is  judged  to  be  similar  to  that  of  convenient  test  courses  (e.g.,  data  from  Aberdeen  Test 
Center  Munson  Gravel  course  can  be  used  as  a  replacement  for  the  Churchville  “C”  course). 

c.  An  ex  posure  t  ime  f  or  each  t  est  s  peed  ( by  t  est  course)  w  ill  be  de  tennined  f  rom  a 
measured  probability  distribution  function  of  a  recent  endurance  test  of  a  similar  vehicle  or  will 
be  computed  using  the  assumption  of  a  beta  distribution.  The  beta  distribution  is  defined  as: 
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Probability  Density  Function 


Xa  ^l-X)13  1 


f(X)=  r  . 

V  '  B[a,|3] 

where  B[a,(3]  is  the  beta  function 
with  parameters  a  and  (3,  given  by 
1 

B[a>(3]  =  Jxa_1(l-X)(3_1dX 


d.  The  parameter  alpha  is  based  on  the  ratio  of  the  average  speed  to  the  maximum  speed 
and  is  selected  from  T able  3.  T  he  parameter  beta  is  calculated  using  an  optimization  routine 
such  t  hat  t  he  r  esultant  av  erage  s  peed  eq  uals  t  he  d  esired  av  erage  s  peed.  An  e  xample  us  ing  a 
minimum  speed  of  10  km/hr  (6  mph),  an  average  speed  of  23  km/hr  (14  mph),  and  a  m  aximum 
speed  of  32  km/hr  (20  mph)  is  shown  in  Figure  6. 


TABLE  3.  VALUES  OF  ALPHA 


RATIO  OF  AVERAGE 

TO  MAXIMUM  SPEED 

ALPHA 

10% 

0.2 

20% 

0.5 

30% 

0.75 

40% 

1.25  I 

50% 

2 

60% 

3 

70% 

4 

80% 

6 

90% 

13 
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Example  Beta  Distribution  of  Speed 
Percent  of  Time  at  Speed 


0  5  10  15  20  25 

Speed  -  mph 


Figure  6.  Example  beta  distribution. 


e.  The  exposure  times  for  each  test  speed  are  combined  with  the  corresponding  ride  level 
to  compute  an  overall  weighted  acceleration  using  (equation  B3  from  ISO  2631-1): 


a 


ws 


1/ 

72 


where: 

aws  =  equivalent  vibration  magnitude  (rms  acceleration  in  m/s  ). 

awi  =  vibration  magnitude  (rms  acceleration  in  m/s2)  for  exposure  duration  T, 

Ti  =  exposure  duration,  minutes. 

f.  The  overall  weighted  acceleration  can  be  used  to  provide  a  health  hazard  assessment 
by  terrain  type  (primary,  secondary,  off-road). 
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APPENDIX  A.  VEHICLE  INFORMATION. 


Vehicle  Model:  _ 

Year  of  Manufacture/Delivery: 

Serial  Number:  _ 

Odometer  Reading:  _ 

Engine  Hours:  _ 

Type  Test  Load:  _ 

Vehicle  Weight  as  Tested: 

Total:  _ 

Front  Axle:  _ 

Intermediate  Axle:  _ 

Rear  Axle:  _ 

Tire  Information: 

Front: 

Manufacturer:  _ 

Size:  _ 

Type:  _ 

Pressure:  _ 

Rear: 

Manufacturer:  _ 

Size:  _ 

Type:  _ 

Pressure:  _ 

Seat  Description: 

Driver:  _ 

Passenger:  _ 

Seat  Setting: 

Driver:  _ 

Passenger:  _ 

Condition  of  the  Vehicle: 
General  Comments: 
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APPENDIX  B.  WBV  TEST  COURSES  AND  OPERATIONAL  MODE  EXAMPLES. 

TABLE  B-l.  TEST  COURSES  AND  SPEEDS  FOR  WBV  ANALYSIS. 


TERRAIN 

TYPE 

COURSE 

APPROXIMATE 
ROUGHNESS 
(in  rms) 

TOP 

SPEED 

(mph) 

SPEED 

INCREMENT" 

(mph) 

Primary 

Perryman  Paved 

0.07 

50 

10 

Secondary 

Perryman  1 

0.3 

35 

5 

Secondary 

Perryman  A 

0.3 

35 

5 

Cross  Country 

Perryman  2 

0.8 

25 

5 

Cross  Country 

Perryman  3 

2.75 

20 

5 

a  Test  Engineer’s  discretion. 


TABLE  B-2.  EXAMPLES  OF  WBV  ANALYSIS  FOR  OPERATIONAL  MODE. 


VEHICLE  TYPE 

VEHICLE  CAPABILITIES 

DATA  NEEDS 

REASONING 

High  Mobility 
Engineer 

Excavator 

Earthmoving  machine. 
Performs  backhoe,  loader, 
and  bull  dozing  operations. 

Collect  data  while 
vehicle  is  performing 
bull  dozing  and 
backhoe  operations. 

Multiple  shock  forces 
from  buckets  engaging 
in  operations  can 
transfer  to  seat. 

XM1157  10-ton 
Dump  Truck 

Transport  and  spread  up  to 
20,000  pounds  of  aggregate. 

Collect  test  course  and 
speed  data  while 
vehicle  is  loaded  and 
unloaded. 

Weight  of  aggregate 
will  affect  WBV 
exposure  due  to  added 
weight  on  suspension 
system  of  vehicle. 

Airborne  Scraper 
and  Water 
Distribution 

System  (ASWDS) 

ASWDS  performs 
earthmoving  road 
construction  and  water 
distribution  functions. 

Collect  movement  data 
while  simulating 
scraping  operations, 
unloaded  and  loaded 
water  operations. 

Scraping  operations 
will  impart 
significantly  more 

WBV  into  vehicle  cab 
than  the  water 
distribution. 

Bituminous 

Material  Paving 
Machine  (BMPM) 

Performs  paving  operations. 

Simulating  paving 
operations  by  loading 
hopper  and  activating 
extended  screed. 

Added  load  and 
extended  screed  may 
shift  center  of  gravity, 
resulting  in  different 
WBV  exposure. 
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APPENDIX  C.  ABSORBED  POWER  CONSTANTS. 


Longitudinal  (X) 

K0  4.3532 

Ki  1.356 

W,  27rfi 

Fi  1.0 

F2  0.219106 

F3  -0.0185309Wi2+ 1 

F4  -0.0006 1 893 Wf  +  0.219106 

Note:  f;  =  Center  frequency  (Hz)  of  ith  spectral  band. 


Transverse  (Y) 

K0  4.353 

Ki  1.356 

Wi  27rf; 

Fi  0.24052 1 24x  1  O'3  W4  -  0.066974483W2  +  1 

F2  0.57384538x10'5W4  -  0.50 1 704 1 3x1  O^W,2  +  0.33092592 

F3  -0.14979958xl0'5Wi6  +  0.0010088882W14-0.10108617Wi2+  1 

F4  -0.1713749xl0'7Wi6  +  0.53 13735  lxlO'4W  4  -  0.01 1096507Wf  +  0.33092592 

Note:  f;  =  Center  frequency  (Hz)  of  ith  spectral  band. 
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APPENDIX  C.  ABSORBED  POWER  CONSTANTS. 

Vertical  (Z) 

K0  4.3537 

Kj  1.356 

Wi  27rf; 

Fj  -0.1 0245x1 0'9Wi6  +  0.17583x10'5W14-0.44601x10'2W12  +  1 

F2  0.12882xl0'7Wj4  -  0.93394x1 0'4Wf  +  0.10543 

F3  -0.454 16xlO'9Wi6  +  0.37667x10'5W4  -  0.56104xl0'2W2  +  1 

F4  -0.21  179x1  0'11  W,6  +  0.51728x10'7W4  -  0. 17947x1  O^W,2  +  0.10543 

Note:  f;  =  Center  frequency  (Hz)  of  ith  spectral  band 
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APPENDIX  D.  ISO  2631-1  WEIGHTING  FACTORS. 


TABLE  D-l.  WEIGHTING  FACTORS. 


ONE-THIRD  OCTAVE 
CENTER  FREQUENCY  (Hz) 

WEIGHTING  FACTORS 

LONGITUDINAL/TRANSVERSE  (X,Y) 

VERTICAL  (Z) 

0.2 

0.243 

0.121 

0.25 

0.365 

0.182 

0.315 

0.530 

0.263 

0.4 

0.713 

0.352 

0.5 

0.853 

0.418 

0.63 

0.944 

0.459 

0.8 

0.992 

0.477 

1.0 

1.011 

0.482 

1.25 

1.008 

0.484 

1.6 

0.968 

0.494 

2.0 

0.890 

0.531 

2.5 

0.776 

0.631 

3.15 

0.642 

0.804 

4.0 

0.512 

0.967 

5.0 

0.409 

1.039 

6.3 

0.323 

1.054 

8.0 

0.253 

1.036 

10.0 

0.212 

0.988 

12.5 

0.161 

0.902 

16.0 

0.125 

0.768 

20.0 

0.100 

0.636 

25.0 

0.080 

0.513 

31.5 

0.0632 

0.405 

40.0 

0.0494 

0.314 

50.0 

0.0388 

0.246 

63.0 

0.0295 

0.186 

80.0 

0.0211 

0.132 

100.0 

0.0141 

0.0887 
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APPENDIX  E.  HEALTH  HAZARD  ASSESSMENT  DATA  FILE  FORMAT. 

E-l.  Header. 

Each  BCR  file  starts  with  a  header  that  includes  three  descriptive  groups  of  infonnation  about 
the  test  parameters  (channel  labels,  engineering  units,  sampling  rate,  etc.).  The  filename 
(maximum  of  12  characters)  appears  on  the  first  line  of  the  header.  A  descriptive  title  can  be 
placed  on  the  first  and  second  line  starting  at  column  17.  Header  infonnation  is  placed  in  blocks 
that  start  with  a  label  that  is  surrounded  by  angle  brackets:  <TEST>,  <SAMPLING>,  and 
<CHANNELS>.  Blank  lines  in  the  header  section  are  ignored. 

E-2.  Data. 


The  data  section  begins  with  the  label  <DATA>,  followed  by  a  line  containing  the  fonnat  of  the 
data  (in  Fortran  style)  and  then  the  numeric  data  itself  organized  in  columns  with  each  column 
representing  one  channel  of  data.  The  data  may  be  delimited  by  commas,  tabs,  or  spaces,  with 
one  row  representing  one  time  sample  of  all  signals.  Each  data  point  can  have  a  maximum  of 
seven  characters.  The  best  way  to  accomplish  this  is  to  limit  the  resolution  of  the  data  to 
thousandths  (0.001). 

E-2.1  The  Jolt  program  is  designed  to  treat  each  triaxial  accelerometer  as  a  single  data  file. 
Therefore,  for  each  measurement  location,  there  needs  to  be  a  separate  BCR  data  file  containing 
the  acceleration  data  for  that  location.  The  Jolt  program  is  also  designed  to  read  input  from  each 
of  the  three  axes  (X,  Y,  Z)  of  the  accelerometer  simultaneously.  Therefore,  the  BCR  data  file 
has  to  contain  three  columns  of  acceleration  data  with  the  order  of  X  being  the  first  column,  Y 
being  the  second,  and  Z  being  the  final  column.  Data  may  only  span  512  characters  across  each 
row  of  data. 

E-2.2.  The  program  requires  the  acceleration  data  be  in  the  units  of  m/s".  For  convenience  of 
the  user,  and  for  those  files  in  units  other  than  m/s2  (such  as  g’s),  the  user  may  specify  a 
conversion  factor  and  have  the  program  convert  the  data  to  m/s2.  Figure  E-l  is  an  example 
illustration  of  the  BCR  file  data  structure.  Note:  You  must  include  keywords  as  they  are  shown 
in  red  (including  colons). 
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APPENDIX  E.  HEALTH  HAZARD  ASSESSMENT  DATA  FILE  FORMAT. 


A  Sample. BCR 


<TEST> 

Description: 

Test  Location: 

Vehicle  Type : 

Seat  Position: 

Test  Terrain: 

Vehicle  Speed: 

<SA11PLIHG> 

Humber  of  Channels :  3 

Samples  per  channel:  25298 
Sampling  rate  (Hz):  416.667 

Signal  duration  (s):  60.7 


<CHAH11ELS> 

- 123456- 

-12345678 

01 

seat 

X 

acceleration 

m/sA2 

XAC 

seat 

02 

seat 

y 

acceleration 

m/sA2 

YAC 

seat 

03 

seat 

Z 

acceleration 

m/sA2 

ZAC 

seat 

<DATA>  ************************************************** 
1  (F7 . 2 )  , 2  (F8 . 2 )  [Obsolete  Fortran  format  ...  It  is  ignored] 

0.05  0.21  0.14 
0.05  0.19  0.13 

- etc .  .  . 


THIS  IS  A  SAMPLE  TITLE 

Human  response  to  shock  and  vibration 


Phase  2  -  Vehicle  Data 

outside  track 

SUV 

Driver 

gravel 

15  mph 


Figure  E- 1 .  BCR  data  structure. 
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APPENDIX  F.  ABBREVIATIONS. 


BCR 

British  Columbia  Research 

COV 

Coefficient  of  Variation 

DOD 

Department  of  Defense 

HDBK 

Handbook 

Hz 

Hertz 

ISO 

International  Standards  Organization 

km/hr 

kilometers  per  hour 

m/s2 

meters  per  second  squared 

mph 

miles  per  hour 

MTVV 

Maximum  Transient  Vibration  Value 

NATO 

North  Atlantic  Treaty  Organization 

OMS/MP 

Operational  Mode  Summary/Mission  Profile 

PSD 

Power  Spectral  Density 

rad/s2 

radians  per  second  squared 

rms 

root  mean  square 

SAE 

Society  of  Automotive  Engineers 

TARADCOM 

Tank- Automotive  Research  and  Development  Command 

TOP 

Test  Operations  Procedure 

VDV 

Vibration  Dose  Value 

WBV 

Whole  Body  Vibration 
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